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A reduced fast component-by-component construction of lattice point 
sets with small weighted star discrepancy 

Ralph Kritzinger* and Helene LairrieR 


Abstract 

The weighted star discrepancy of point sets appears in the weighted Koksma-Hlawka inequality 
and thus is a measure for the quality of point sets with respect to their performance in quasi-Monte 
Carlo algorithms. A special choice of point sets are lattice point sets whose generating vector can 
be obtained one component at a time such that the resulting lattice point set has a small weighted 
star discrepancy. 

In this paper we consider a reduced fast component-by-component algorithm which significantly 
reduces the construction cost for such generating vectors provided that the weights decrease fast 
enough. 
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1 Introduction 


Given an IV-element multiset of points {xo,... , x^v-i} £ [0,1), we may approximate integrals over the 
s-dimensional unit cube by a quasi-Monte Carlo (QMC) rule, i.e., 
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N 


N—l 


J2 A x d- 

n=0 


For detailed information on QMC-integration see mmmm- 

In 1998 Sloan and Wozniakowski m introduced the concept of weighted function spaces where 
each group of coordinates is equipped with some weight according to its importance. Denote the set 
{ 1 ,..., s} by [s] and let 7 = (7 u ) u c[s] be a weight sequence of non-negative real numbers, which model 
the importance of the projection of the integrands / in the weighted function space onto the variables 
Xj for j € u. A small weight 7 U means that the projection onto the variables in u contributes little to 
the integration problem. In the present work we consider a special choice of weights, so-called product 
weights (~fj)j >\, where 7 U = TljeuTi an d 70 := 1 , and in particular, the weight 7 j is associated with 
the variable xj. 

In this paper we assume that 7 = ( 7 j)j>i is a non-increasing sequence of positive weights with 
7 j < 1 and ( 7 u )uc[s] are the corresponding product weights. Such weights are useful when considering 
functions whose dependence on successive variables is decreasing. 


A particularly important kind of point sets for QMC-integration are so-called lattice point sets. 
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Research Program “Quasi-Monte Carlo Methods: Theory and Applications”. 
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They originated independently from Hlawka [S] and Korobov m- a lattice point set H/y(z) — 
{xo,..., Xjv_i} can be constructed with the aid of a generating vector z. For a positive integer N > 2 
and a vector z £ {1,..., N — 1} S the corresponding lattice point set is of the form 

Pat(z) ={{A z ] : k= lj . 

The brackets {.} around -^z indicate that we take the fractional part of each point. For vectors, {.} 
is applied component-wise. See m EU22]. 

We want to measure the quality of lattice point sets Pn( z ) with respect to their performance in a 
QMC rule. Therefore we define the weighted star discrepancy. 

Definition 1. Let 7 = ( 7 u ) u c[sl be a weight sequence and let Pjy = {xo,..., xjv-i} C [0, l] s be an 
N-element point set. The local discrepancy of the point set Pjv at x = (aq,... ,x s ) £ [0, l] s is defined 
as 

1 s 

discr(x,P/v(z)) := — ^ X[o, x )(p) ~ II X L 
pePjv 3 =1 

where X[o.x) denotes the characteristic function of [ 0,x). The weighted star discrepancy of P^ is then 
defined as 

D* n JP N ) ■■= sup max 7 u |discr((x u , 1 ), P/v)|. 

x€(0,l] s 

We denote the weighted star discrepancy of a lattice point set corresponding to some generating 
vector z by D* Nl (z), as this Pn{ z) is completely determined by z. To see why the weighted star 
discrepancy is a measure for the quality of our point sets we study the following identity of Hlawka 
[7] and Zaremba gl] (see also HIED, given by 

QnM)~ i sU)= J2 (“ 1 ) |U| 7 U / discr((x u , 1 ), -P/v(z))7 u vr / ( x u> 1 ) dx u , 

0 ^uc [s] 7 M |u| dx u 

where QN,s{f ) = Sj=i /( x j) denotes the QMC-rule, I s = Jj 0 X j S /(x) dx the integral operator and 
(x u , 1 ) the vector (aq,..., x s ) with Xj = Xj if j £ u and Xj = 1 if j f u. 

Applying Holder’s inequality as in HE2] for integrals and sums we obtain 

|Qiv, s (/)-/ s (/)|<^, 7 ( Z )||/|| 7 , (i.i) 

where ||.|| 7 is some norm dependent on 7 but independent of the point set Pjy(z). If / is sufficiently 
smooth ||/|| 7 coincides with the weighted variation of / in the sense of Hardy and Krause. The first 
factor in ( 11 . 11 ) is the weighted star discrepancy of the point set P/v(z) and depends only on P/v(z) and 
the weights. Thus we see that the quality of a lattice point set Pn{ z ) is the better the smaller its 
weighted star discrepancy D* Nl (z). We want to find lattice point sets Pn( z ) with small weighted star 
discrepancy. 

As no explicit constructions for good lattice point sets are known for dimensions s > 2 , one usually 
employs computer search algorithms to find good generating vectors. There exist many papers on the 
construction of generating vectors for lattice point sets with a small weighted star discrepancy: Joe 
[9] has given a component-by-conrponent construction for generating vectors of lattice point sets with 
a prime number N of points, which have a weighted star discrepancy of order for any 6 > 0. 

Their generating vector has a construction cost of order sN log N, where an approach of Nuyens and 
Cools m can be used to reduce the construction cost. 

In [2U] Joe and Sinescu have achieved the same results for a composite number of lattice points 
and product weights. Finally in m they considered general weights and a prime number of points. 
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Dick et al. [2; have given a reduced fast algorithm for the construction of generating vectors of 
lattice point sets with N a prime power. They varied the size of the search space for each coordinate 
according to its importance and considered the worst-case error of integration in a Korobov space to 
measure the quality of their lattice point sets. 

Let b be an arbitrary prime number and m a positive integer. In the present work we consider lattice 
point sets with N = b m elements and study their weighted star discrepancy. As mentioned before, the 
generating vector z = (z \,..., z s ) of such lattice point sets can be obtained one component at a time. 
When using the standard component-by-component construction, in the following frequently abbrevi¬ 
ated by CBC construction, each component is chosen from {z G {1,2,..., b m — 1} : gcd (z, b m ) = 1}. 
As done in [2] for the worst-case error, we speed up the construction of such generating vectors by 
reducing the search space for each component, while still achieving a small weighted star discrepancy 
of the corresponding lattice rule. To this end we define non-decreasing 0 < w\ < W 2 < ... G N and 
set 
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N,Wj 


{z G {1,2,... , b m Wj — 1} : gcd (z, b m ) = 1} if Wj < m, 
{1} if Wj > m. 


Note that these sets have cardinality b m ~ Wj ~ 1 (b — 1), for Wj < m. In what follows we denote by Zfj w 
the cartesian product b wl ZN, Wl X ... X b Ws Z]y tWs , where b w iZ^^ Wj means that every element of Z\r Wj is 
multiplied by b w 'i. We denote by z € Z S N w a vector z = (b wi z\, ..., b Ws z s ), with Zj G Z^ Wj for j G [s]. 
We study the weighted star discrepancy of lattice point sets P/\r(z) with generating vectors z G Z' S N w . 
Dick et al. J2l have considered the worst-case error for approximating the integral of functions in 
suitable spaces by a QMC rule based on lattice point sets. Here, in contrast, we study the weighted 
star discrepancy of these lattice point sets which is another important quality measure. We will see 
that for sufficiently fast decreasing weights we can construct lattice point sets with small weighted 
star discrepancy, while significantly reducing the construction cost in comparison to the standard CBC 
construction. 


It follows from m Theorem 3.10 and Theorem 5.6] that 

D* n ^{z) < W U ~ f 1 - +^P^ l7 (z), 

uC[s] V V 

where 


( 1 . 2 ) 


= 51 7uPtv(z,u) 
uC[s] 


(1.3) 


and 


N -1 


/ 


r N (z,u )=—n 


k =0 j £u 


E 


V 




^2iT\hkb w j Zj /N 

\h\ 


- 1. 


0 


(1.4) 


Using this estimate for the weighted star discrepancy we derive the results in Sections [21 El and 01 
Finally, we introduce the concept of tractability [15, 16, 17 . To this end we define the information 
complexity (often refered to as inverse of the weighted star discrepancy) as 


N*(e, s ) = min {A" G N 0 : D* N ^(z) < e}, 

which means that N*(e,s) is the minimal number of points required to achieve a weighted star dis¬ 
crepancy of at most e. Of course we want the information complexitiy to be as small as possible. 
Therefore we are interested in how fast it increases when e^ 1 and s grow. We define the following 
notions of tractability. We speak of 
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• polynomial tractability, if there exist constants C. t\ >0 and 72 > 0 such that 
N*(e, s ) < Ce~ T1 s T 2 for all e G (0,1) and all s G N and of 


• strong polynomial tractability, if there exist positive constants C, r such that 

n(e, s) < Ce~ T for all e G (0,1) and all s € N. 


Roughly speaking, a problem is considered tractable if its information complexity’s dependence on 
e _1 and s is not exponential. We will show that the above mentioned reduced fast component- 
by-component construction finds a generating vector z of a lattice point set that achieves strong 
polynomial tractability if 

OO 

J2^ bWj < 00 

3 = 1 


with a construction cost of 


min{.s,t} 

O ( IVlog N + min{s,t}N + N ^ (m — Wd)b 

d =1 


-w d 


operations, where t = max {j G N : Wj < m}. 

The structure of this paper is as follows. In the next section we derive an upper bound for the 
arithmetic mean of the weighted star discrepancy over all possible lattice point sets constructed by 
a generating vector z G Z s Nvj . In Sections [3] and U we present a reduced fast CBC construction for 
generating vectors of lattice point sets with small weighted star discrepancy. Finally, in Section Owe 
study conditions on the weights 7 j and wj for achieving strong polynomial tractability. 


2 The arithmetic mean over all z £ Z s Nyv 


First of all we estimate the arithmetic mean of the weighted star discrepancy over all possible gen¬ 
erating vectors z = (b Wl zi ,..., b Ws z s ) G Z S N w , proceeding similarly to [2] and [20j- This yields the 
existence of a lattice point set with small weighted star discrepancy. The upper bound which we ob¬ 
tain for the arithmetic mean is not the same as for the reduced CBC construction in the next section. 
Nonetheless, we need large parts of the calculation of the present section to obtain the estimate in 
Section 0 


Theorem 2.1. Let N = b m , {wj)j> 1 and G Z S N w be as above and let m > 5. Then there exists a 
generating vector z = (b Wl zi ,..., b Ws z s ) G Z S N ^ whose corresponding lattice rule has weighted star 
discrepancy 


D* N)1 { Z)< £ 7u(l- fl¬ 


ue S 


l \ \ M l / l 


N 


+ 


2 N 


UWi + ySN) 


3 =1 


m— 1 


n (Pi+nSN) n *-nft]. 


p=0 


3 = 1 

Wj >m—p 


3=1 

Wj<m—p 


3=1 


with /3j = 1 + 7 j for all j G N and 


5 -= £ w 


( 2 . 1 ) 


h^O 


4 


Remark 1. Provided that the q j’s are summable the bound in Theorem A2.1\ is of order IV 5 log IV for 
arbitrary 6 € (0,1) with an implied constant independent of N and s. Furthermore note that if all 
weights Wj = 0 then we obtain the result in m Theorem 1 and Corollary If. 


Proof. As the first sum in (1 1.2 f) is independent of z, it is obviously enough to consider the mean 
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Mn,s ,7 | ?s | 55 Rn,-/{ z ) 


J N,w\ Z& z s . 


( 2 . 2 ) 


of the second sum. 

We have from 0 p. 186, eq. 9] 


N—l s 






k =0 j =1 


p j + 7 i E 

-<h : 
h^O 


-f< h <f 


^2 nihkb w 3 Zj /N 

\h\ 


- n ft 

3 = 1 


1 s -j N—l s 

= F n(ft+vft) + -x;n 


3 = 1 


JV 


k =1 1=1 


P3 + 7? 

A/0 


^2nihkb w 3 Zj /N 


- n ft- 

l=i 


Thus 


s AT—1 s 

=77 n (ft+vftv)+77 e n 


(2.3) 


5=1 


N 


k=l j=1 


'N,w 
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3 Zj £zZ]] W j 


/ 


\\ 

Pi + 7 i 

E 

^2'Kihkb 3 Zj /N 

\h\ 


-4</x<4 

V 

h^O 
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-n ft 

3 = 1 


s -j AT—1 5 

l=i fe=i l=i 
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l<h< 
"h^ 0 


iV ’^' Zj£ Z N ,Wj -%<h <f 


^27 rihkb w j Zj /N 


To avoid lengthy formulas we use the following abbreviations: 

p2nihkb w 3 Zj /N 


T N , Wj (k) := 55 55 

}<h : 


Zj£Z N , Wj -f<h<% 


\h\ 


and 


i N—l s / \ 

L ».‘.y ■= 17 E II ft + tP—.TnmW . 

iV fc=l j=l V \^N, Wj \ J 


fc=l 3=1 

Then we have 

2 s s 

H (Pi + li^N) + Ln,s,~i — UPr 

l=i l=i 

We study T]y tWj (k) distinguishing the two cases Wj > m and Wj < m. 


- n ft- 

1=1 


(2.4) 


(2.5) 

( 2 . 6 ) 


Case 1: wj > m. This yields Z^ tWj = {1} and thus 

J2'K\hkb w 3 /N 

TN, Wj {k) = 55 j^i = E 


Jlirihkb i 


WA—m 


0 


h^O 


\h\ 


~ T Wr ^ N . (2.7) 


h^O 
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Case 2: Wj < m. Then Z^ tWj = {z £ {1,2,..., b m ~ Wj — 1} : gcd (z, N) = 1}. According to fl2.5li we 
have to calculate Tpf >Wj (k) only for k £ {1,..., b m — 1}. We display these k as k = qb m ~ w i + r with 
q £ {0,..., b Wj — 1}, r 6 {0,..., b m ~ Wj — 1} and (q, r ) / (0,0). Then 


Tn,Wj (k) 
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1 

W\ 


E 

Zj 


^27ri h^qb 171 w o -\-r)b w i Zj /N 


E 


1 

W\ 


E 

£%N,Wj 


^Znihrzj/b m w i 


E 


1 

N 


E 

Zj(z:ZN,Wj 


^2 irihqzj ^Znihrzj /6 m 


( 2 . 8 ) 


If r = 0, i.e. fc a multiple of Wj , this yields 


TN, Wj {k) = E TTT E 1 = |-2W,^|<Sjv- 

~i<h< f 1 1 z j£^N,Wj 

0 


(2.9) 


Next we investigate r £ {1,..., b m Wj — 1}. For any z 3 £ {0,..., b m Wj — 1} we find gcd (zj, N ) = 
gcd (zj,b m ~ w i) £ {1, b, b 2 ,..., and hence 


E k-(d) = E Md) 

d\ gcd ( z j , N ) d| gcd ( Zj ,b m ~ w i ) 


1 iff gcd (z.,-, IV) = gcd (zj, b m Wj ) = 1, 
0 otherwise, 


where /x denotes the Mobius function. 

For any z 3 £ {1,..., b m ~ w 3 — 1} this implies Zj £ Z^ yWj if and only if ]T ^id) = 1. 

d\ gcd (zj ,b m ~ w j ) 

Inserting this fact into (|2.8I) we have 


Tn, W j (k) 


E 

jv w,. 


1 

]h\ 


jrn-Wj _j 

^2-Kihrzj /b m ~ w i 

Zj = l 


E did)- 

d\ gcd (zj,b m ~ w i) 


( 2 . 10 ) 


Studying the two inner sums we find 


» ni — uj a -| 

0 d —1 


£ e 2.i hrzi/b—i £ ^ d) 

z i ml d\ gcd [z j ,b m ~ w i) 


j^rn—w 


j -l 


£ //(d) £ e 2 nihrzjlb—i 




«j=l 

d|2,- 


£ //(d) E e 2M / b ^ W \ 

d\b m - u 'i a=1 


r im — Wj 'j _ 

where the latter equality holds since a G < 1,..., —^1 yields ad G {d, 2d,..., b m ~ w -'} 
= {1 < z,- < - 1 : d\zj} U {b m ~ w i} and 


E d( (l ) = 0 , 

d\b m - w j 


( 2 . 11 ) 


since w 3 < m. 
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Changing the order of summation we obtain with ( 12 .lip 


l ~ w j —i 

E e— 


2'K\hrzj /6 m w i 


E M= E 4^)£e 2 * ihra/d = E 


d|gcd(2:j,6 m w 3 ) 


d\b m ~ w j 


a= 1 


d\b rn ~ w o 

d\hr 


frm-Wj 

d 


With (I2.1UP this leads to 


1 ( h m ~ w i 

TN, Wj (k) = E |^| E d ^\ 


-f <h<% ' 1 d\b m - w 3 
0 d\hr 

d 


d 


= E 

d\b m ~ w j 


' b m- Wj \ 1 

~ZT~ I E Jh\- 


-f<h<f 

h^O 

d\hr 


Using that d\hr is equivalent to gcd ?^ r ) \h we display TN, Wj (k) as 


Tn,Wj (k) = E d f JL \ 

d\b m - W 3 


b m-w j 

d 


E 


gcd (d,r) 


w 


( 2 . 12 ) 


For further investigation of Tjy jWj (k) we first study sums of the same type as the inner sum in 
f!2.12p . For any positive integer a we have 


„£„isr 


(2.13) 


/i^O 

a|/i 


where SV is defined analogously to m- Combining (12.131) with ( 12 . 121 ) we obtain 


T N , Wj (k) = E dn 


b m w .i \ gcd (d, r) 


tfn-vjj 


5 fgcd (d,r)= E ^ g cd ( d ^) 5 fgcd(d,r) 


d 


d 


-gcd (b m gcd (prn-Wj r ) gcd (b m ^ V)Vi+i gcd (r—”i-i,r) 

= b (^S^wj+u ^+1), 


(2.14) 


with € { 0 ,... , m — Wj — 1 }. 

Summarizing, we have for k G {1 , ,b m — 1} 


TN, Wj {k ) = < 


if Wj > m, 

if Wj < m and k = 0 (mod b m ~ w 3 ), 


Sn 

\Z N , Wj \S N 

b (^S^wj+u 5)../!,++1) 

with b v = gcd ( b m ~ w 3 , r) if Wj < m and k ^ 0 (mod b m ~ w 3) 


(2.15) 


Let us choose t G No such that Wj < m for all j < t and wt+i > m. (If t = 0, then Wj > m for 
all j G N. In that case we obtain the generating vector z = (b Wl , ..., U" s ).) With this we are able to 
write Ln , s ,7 from formula ( 12 .5p as 


N—l min{t,s} / \ s / > 


1 


fe = l j=l 


- / N,Wj 
N—l min{£,s} 


j=t+l 


■'N.Wj 


(2.16) 


- t? n (#?■+e n (+ 

i=i+i fc=i j=i V 


7i 


7 T N)W Ak) . 
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Next we aim at finding bounds for '\ f° r w j < m - 

If k is a multiple of b m ~ wj we see immediately from (12.15 p that 

Tn,Wj ( 6 ‘) .w.j | S N 

—r = ~r?—n = bN ■ 

| N,Wj | | ^- , N,Wj | 

If k is not a multiple of b m ~ w 1, we use a formula from Niederreiter m for S n with arbitrary n £ N, 
given by 

S n = 21ogn + 2 y - log 4 + e(ra), (2.17) 

( l | \ 

where 7 denotes the Euler-Mascheroni constant 7 = lim X] t — log/ ~ 0.577216 ... and 

l -+°° \k =1 / 

j — -j < e(n) < 0 , if n is even, 


“4 < e{n) < 


if n is odd. 


(2.18) 


From (12.151) we know 


TN,wj(k) — b v (S,wj+v — S,wj+v+i) < 0. 


(2.19) 


With m > 5 we find —2 < — y < 0 for wj < m and k not a multiple of b m Wj as follows. The 

upper bound follows immediately from (12.191) . It remains to show the lower bound. First we consider 
Tn.Wj (k) using (12.171) . We have 

T N}Wj (k) = b v {S hV>j+v - S bWj+ , + 1 ) = 6" (—2log b + £(b Wj+u ) — e(b Wj+u+1 )^ 

= -2b v \ogb + b v (e(b w i +v ) - e(b w i +u+1 )) . 

With (12.181) we obtain 

b v (e(b w i +u ) - e(b w i +v+1 j) I < I b v (e(6^ +I/ )) I + I b v (e(b w * +u+1 j) 


< 4 b~ 2u ’i- u ( 1 + 


6 2 


Thus 


T N ,wAk) ^ b Wj-m+l 


> -- 


■>N,Wj 


b- 1 


-26 v log6- 


b- 1 


4 b~ 2w i~ u ( 1 + 


b 2 


Recall from (12.151) that v = log fe (gcd ( b m Wj , r)) £ {0, 1,... ,m — Wj — 1}. Thus 


Tn,Wj (fe) ^ 2^ Wj—m+l+m—Wj—l ^ ^ 




6-1 

Now, with the assumption m >5, 


6-1 

1 
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1 + 
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1 + ^ 


T N,Wj(k) > _ 2 ]og6 _ 4 ^_5 _|_i_1 
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1 

1+ 62 


> -2 lc> " 2 - 4 • 2 _5+1 
2-1 


1 + 


> - 2 , 
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and hence 


-2 < 


Tn,Wj (k) 

| - 2 N,Wj | 


< 0 for wt < m and b n 


t k. 


For any integer p £ {0, ... ,m — 1} with bP \ k and bP +1 \ k the condition b m Wj \ k is equivalent to 
m — Wj > p or Wj < m — p, respectively. Thus we can display (12.161) as 


1 

Ln,s,7 = T7 II (Pj +7 jSn) 

j=t +1 

m—1 N—l min{t,s} / 

x e e n U+ 

p =0 k=l j =1 \ 

Wj>m—p 

bP +1 \k 
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min{f,s} 


r^W*) n & + 


7j 


3 = 1 

Wj<m—p 


\ ^N,Wj 


-\T N)W Ak) 
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JV—1 min{t,s} 


min{£,s} 


n (ft+T e n w+v-sa,) n a-. 


j-t+i 


p=0 fc=l j=l 

i^lfc Wj>m—p 
bP+ 1 \k 


3 = 1 

Wj <m—p 


Tn,wa (&) 


where the latter estimate holds since /3j > 1 , —2 < 3 | < 0 and 7 j < 1 


. From 


|{fe€ {1,..., IV — 1} : 

= |{fe G { 1 ,..., b m - 1 } : W | k}\ - |{fc £ { 1 , • • •, b m - 1 } : lf +1 \ k} 
= b m ~ p - 1 - ( 'b m ~ p - 1 - l) = b m ~ p ~ l {b - 1 ) 


we get 


m—1 


min{£,s} 


min{t,s} 


^,7<^ n ^-+7 A)e^-i) n (P 3 + 13 SN) n &■ 


j=t +1 


p =0 


i=i 

Wj>m—p 


3 = 1 

Wj<m—p 


Inserting this into (12.61) we obtain for the arithmetic mean 
1 


Mn,s ,7 — jy IT (^3 7 j'S'iv) 


i=i 


m—1 


min{i,s} 


min{£,s} 


( 2 . 20 ) 


n (ft-+7^) n &-n&- 


( 2 . 21 ) 


j=t+i 


p=0 


i=i 

Wj>m—p 


3 = 1 1 =1 

Wj <m—p 


This proves with ( 11 . 21 ) the existence of a vector z £ i?^r w such that the weighted star discrepancy 
D * N)1 { z ) fulfils 
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( 2 . 22 ) 
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3 The reduced CBC construction 


In this section we give a component-by-component construction for the generating vector and an upper 
bound for the weighted star discrepancy of the corresponding lattice rule. 

Algorithm 1. Let N = b m and (wj)j >i be as above and construct z = ( b Wl z \,..., b Ws z s ) G Z S N w as 
follows: 

1. Set z\ = 1. 

2. For d G [s — 1] assume z±,... ,z d to be already found. Choose z d+ \ G Zn jW such that 

R d +^z 1 ,...,b w ^z d ,b w ^z) 

is minimized as a function of z. 

3. Increase d by 1 and repeat the second step until z = (b Wl zi ,..., b Ws z s ) is found. 

In the algorithm above the search space is reduced for each coordinate of z according to its 
importance. This results in a considerable reduction of the construction cost as we will see in Section 
U This is why we call this algorithm a reduced CBC-algorithm. 

The following theorem gives an upper bound for the figure of merit Rfj of lattice point sets with 
generating vectors obtained from the algorithm above. 

Theorem 3.1. Let z = (b wi z±,... ,b Ws z s ) be constructed according to Algorithm QJ Then for every 

d € [s], 

< 7 (^i ,..., b Wd z d ) < 1 n (Pj + (l + 2& min{ ^’ m >) 7 jS N ) • (3.1) 

i=i 

Corollary 3.2. Let N = b m and {wj)j>\ be as above and let z = (b wi zi ,..., b Wa z s ) G Z S N w be 
constructed with Algorithm Q] Then the corresponding lattice rule has a weighted star discrepancy 

SJvrW < £ 7» (i - (l - ^)) M + ^ ft (ft + (l + 26"'"7iftv). 

uC[s] V V JV// j =1 

Proof. Combining (11.21) . (12.11) and Theorem 13.II we immediately obtain the result. □ 


To prove Theorem 13.11 we use the the following 

Lemma 3.3. Let N = b m , (wj)j >i and Z^^ Wj be defined as above and recall from (12.41) the notation 

Jl , K\hkb w i Zj /N 

T N, Wj (k) = E E ITT • 

Zj eZ NtW . <h <f 11 

h^O 


Then 


N -1 


E 


\T N , Wj (k)\ 

I -2- N,Wj | 


2^min {uij ,m 


^Sn- 


(3.2) 
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Proof. As before, we distinguish the two cases Wj > m and Wj < m. 
Case 1: Wj > m. Then (|2.15|) yields 


wW)I 

h \ Zn ^\ 


N -1 

J2 S N = (N -1)S N < 2NS n = 2 b min{w i’ m} S N . 

k =1 


Case 2: Wj < m. We use (12.15 D and (12.81) to find 


v^ 1 \T N , Wj (k)\ , ^ 

E \ 7 | + E 

fc=i k =i 

For any r G {1,..., b m ~ Wj — 1} the condition gcd (r, b m ~ Wj ) = 6" is equivalent to b v | r and b u+1 \ r 
simultaneously. Using this we investigate the last sum in the above equation: 


yMI 

h \ Zn ^\ 


\TN, W j(k)\ 
P'N.v: :j | 


w j -1 


= (b w * -l)S N + b w i Y, 


\T N ,wAr)\ 


r =1 


- / N,Wn 


ym—w 


3 


-1 


E 


\TN, Wj (r)\ 

\^N,Wj | 


1 

^N,Wj 


m—Wj — 1 w j — i 

E E PWr)|. 

v=0 r =1 

&"|r 

&"+V 


Once more with the aid of (12.151) this yields 


m—Wj — 1 w j — i 

E E ^(Sywj+v — Sywj+v+l) 


itti—wa -t 

b \T N , Wj (r)\ 1 
r=l l^^'l ~o 7=1 

b v \r 
b u+1 \r 

m—Wj — 1 j/ rL ~ w j _i 

E E ^(V^+i 


1 


12, 


v= 0 r=l 
6"|r 

b"+V 


- S.Wj+v). 


Analogously to (12.201) we find 

{r € {1,..., b m ~ w i - 1} : b u | r and 6 y+1 f r} 


= b 


m—Wn—u—1 


(6-D 


and hence 



r=l 


( r ) I 

\^N,Wj | 


m—Wj — 1 

^=0 


Altogether we have 


TV-l 


E 


\^N,Wj | 


(6 w i - 1)5 jv + b w *(S N - S b »i) < 2 b w *S N = 2b min{w >’ m} S N 


and the proof is complete. 


□ 


With the aid of Lemma 13.31 we are able to prove Theorem 13.11 using induction on d. 
Proof. According to Algorithm [T] we set z\ = 1 in Step 1. We have to show that 

RN^b Wl ) < ^ (ft + (l + 26 min l™ 1,m l) 7l S N ) . 
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With (12.31) we have 


RjV, 7 


(, b wi ) 


1 

N 


N -1 


E 


/ 


0i + 7i E 

~T< h <f 

h^O 


e 2nihkb w i /N 

N 

/ 


-0i 


i 

iv 


1V-1 


E 71 


E 


Q 2-nihkb w i /N 

\h\ 


Again, we consider the two cases wi>m and w\ < m separately. 
Case 1: w\ > m. Then 


Rn, 7 


(b wi ) 


N -1 


M E71 E 


Jlirihkb™-*--™ 


N 

1 

iv 


fc=0 -f<ft<f 
h^O 


\h\ 


= -^JiNS n < y(i + 7i + 2N^iSn) 


(0i + 2& mMm ’™} 7l £ 7V ) < ^ (0i + (l + 26 min ^’ 1 ’ m >) 7 i 5jv) 


which is the desired result. 

Case 2: w\ < m. After interchanging the two sums, once more, we split up the inner sum as follows: 


< 7 (^) 


- E 

N ^ 

ft^O 


1 

W\ 


N—l 

e 27ri/ifc /b Tn ~ w ^ 

k =0 


— V — V e 

jV ^ l/il ^ 
-f <ft<f 1 1 fc=o 

ft^O 

f ) m-w 1 | ^ 


1V-1 


2'Kihk/b rl 


N -1 


11 + — y^ — y^ e 


2 , K\hk/b‘ n 


~f </i<f fc =° 
ft^O 

b m -“ifft 


Now we are able to compute the inner sums. The first one sums to A r , whereas the second one equals 
zero which can immediately be seen by applying the formula for finite geometric series. Thus 


R 1 N , 1 {b w ') = n £ 

0 

b m ~ w 11 ft 


1 

IE 


We use (12.131) to find 
?i 


< 1 (ft + 26 « 7iSn) 

1 (0i + (l + 26 min ^- m >) , 


< 


N 


as claimed. 

Let d £ [s — 1] and assume that we have some z £ such that 

R*r„(b Wl zi, • • •, b”*z d ) < E fj (ft + (l + 26 min (---}) 7 ^iv)• 

l=i 

We have to prove the existence of a z d + 1 £ -2ftv,«j d+1 with 

-| C?+l 

i$g(& wi *i, • • •, b w ^z d+ 1 ) < - n (01 + (l + 26 min ^’ m >) 7j S N ). 

l=i 
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Using again (12.31) we have for any z d + 1 G Z^ )Wd+1 that 
R d ^{b Wl zi ,..., b Wd+1 z d+ 1 ) 


JV-l d 


/ 


n ^ n 


fc=o i=i 


Pj + r/j E 


V 

= Pd+iR d N ,- y (b wl z 1 ,...,b Wd z d ) 


~T< h <T 

0 


^2'K\hkb w o Zj /N 

i*r 


\ 


/ 


Pd+1 + 7d+l X! 

- <h ; 
/i^O 




=l 27ri/i/c6 u;£i + 1 ,2^+1/AT 

N 


- Am n a* 

j=i 


+ 


7d+l 

N 


N-l d 

En 

k =0 j=l 


Pj + 7j E 
~T< h <f 

h^O 


^2TT\hkb w 3 Zj /N 


\ 


\h\ 


E 


I ~f< h <T 


^2nihkb Wd + 1 Zd+i/N 

\h\ 


= f3 d+1 R d NtJ (b w ' Zl ,..., V°* Zi ) + n (ft + 7^iv) 


i=i 


+ 


7d+l 


JV-1 


1 27ri/i.fc6 u ’ d + 1 z ( j_|_i/A r d 

N Ei -f<h<f ^ 3 

h^O 


( 


E E 


A? + 7j E 

-<h-_ 
hjt 0 


-T< h <T 


^2nihkb w 3 Zj /N 


(3.3) 


Next we consider the arithmetic mean of R d Nl (b Wl z \,..., b Wd z d , b Wd+1 z ) over all 2 G Zn, w<1+1 - As 
only the third summand in (13.31) depends on the (d + l)-st coordinate it suffices to investigate the 
mean of this summand. Clearly, if we have some upper bound for the mean over all 2 G Zj^ tWd , there 
exists a z d+ 1 G Z^, Wd+1 which satisfies this bound. Thus we study 


1 


E ^E E 


N—l ^nihkb™ d + 1 z /N d 

\ z€Z~ Wd+1 N mi ^ fl 

0 


( 


pj + 7 3 E 

-<h: 
h^O 


-T< h <T 


-j27 nhkb w o Zj /N 

W\ 


\ 


We bound this term by its absolute value 


1 


-'N,w d+1 1 zGZn 


E 


7d+l 

N 


"d+i 


Af—1 ~27rihkb w d +! z / N d 

- -TT 

* w M 


( 


E E 

-<h< 
h^O 


k=i-f <h <^ 


Pj + 7 ■ 3 E 


V 


A/0 


3 2 7rihkb w i Zj /N 

W 


JV-l T 
< 7d+i 1 


^ fc =1 I^N^d+il 


E E 

^JV lU , d+1 -f <h<f 


=1 27ri/ifc6 u;d + 1 2:^+1 /AT 


N 


d 

><n 

i=i 


( 


Pj + 7? E 


-^2nihkb w 3 Zj /N 


V 




W 


7d+l |Av,i/j d+ i(^)l 77 , c s 

~tr E —r 11 (ft + 7i^) 


< 


^ fc=l I^N^d+il J=i 
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d 

< mi 2 n (/?,+ 7 jS N ), 

3= 1 

where the last estimate stems from application of Lemma 13.31 Combining this with (13.31) we have 
shown the existence of a z d + i € Ziy jWd+1 such that 

R d ^(b w 'z 1 ,...,b Wd z dl b Wd +'z d+1 ) 

< p d+1 R d N ^ Zl ,..., r««z d ) + n (Pi + 7j5jv) + f[ (ft + Tf-W 

j =i i=i 

We use the induction hypothesis to find 

< ^±1 fl (/?,- + (l + 26 min ^’ m >) 7 jS N ) + 7d ^ 5jv n (Pi + 7i-5iv) (l + 26 min ^+!’ m >) 

i=i i=i 

< (flu-i + (l + 2b min ^> m >) 7 rf+ i57v) ^ n (ft + (i + 26 min{ ^’ m} ) 7jS'jv) 

3 = 1 

-i <i+l 

= S n(ft+( i + 2 i > mln{ ” , ' m, )v s ») 

i=i 

which completes the proof. □ 


4 The reduced fast CBC construction 


By now we have seen how we can construct a generating vector of a lattice point set with low weighted 
star discrepancy with a reduced CBC construction as in the previous section. Now we study the 
construction cost of this algorithm. In fact the CBC algorithm can be made faster to construct 
generating vectors for relatively large N and s. To show this we follow closely |2] and m- 

Let d G [s — 1] and assume that we have already found ( b wi z \,..., b Wd z d ). Then we have (cf. (12.31) 1 


N-l d 






k =0 j =1 


Pi + 7 i Y 


V 


h/0 


^lTT\hkb w o Zj /N 

w 


- n ft- 

3 = 1 


Define r(L) = max{l, |/i|}. Then 


Pi + 7? 

/i^o 


3 ! 2nihkb w j Zj /N 


\h\ 


= Pj + 7j Y 


i\ ^ l _IV 


= 1 + 7 i Y 


-f < h <f 


^27ri hkb w o Zj /N 
r{h) 

2iz\hkb w j Zj /N 
r(h ) 


- 1 
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Hence we have 


**,(.... b"’-z d ) = i E n I 1 + H £ 

fc=0 j=l 
TV—1 


TV-1 d 




^2irihkb w i Zj /TV 

r{h) 


-rift 

j=i 


(4.1) 


1 


= jyE - II 


/c=0 


J=1 


where we have defined 


and 


%(&) = II ( 1 + 7/0 

i=i V 

0(®) = 5Z 

~Y< h <i 


kb w iz. 


N 


Jlitihx 


r(h ) 


However, this is exactly the situation as dealt with in m Section 4.2]. Thus we know that 

/ kb w j z ■\ 

(j) ( —at^-) takes on at most N different values, namely 


0(0), 


A 


N- 1 
iV 


which can be computed in 0(N log N) operations and stored in a memory space of size 0(N ), as 
demonstrated in H2j. 

Next we investigate one actual step of the CBC construction. Assuming that we have already 
found ( b Wl zi ,..., b Wd Zd) G Z^ w we have to minimize Rff~^(b Wl zi, ..., b Wd Zd, b Wd+1 z ) as a function of 
z G Zjy jWd+1 to find Zd +1 G -ZTV,io d+1 - For urd+i > m we just set z^+i = 1 and are done. Therefore let 
u>d+ i < m. Considering 631) we have 


1 


TV—1 


d+1 


R d N,^(b Wl zi, • • ■ , &™ d+1 Z d+ i) = — %+i(&) - n Pj 

v fc=o i=i 

1 

' iV 


x>m (i+w(^ ? ^^ ±i ))-nft 

fc=0 V V y y j=l 

4E,^( 1+ w({^}))-nft 

fc=0 V V k } ' ' 7 = 1 


It is obviously enough to minimize jj Vditycj) ({—— \ Zd+1 })■ To do this we proceed analogously 

to [2], We define the matrix 


A = 


fkb Wd+1 z d +i 1 


V 


N 




^d+lG-Z/v,™ • ,fce{0,...,TV-l} 


and 

and find that 


Vd = (Vd(0),r]d{l),---,Vd(N -1)) 
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We can display the matrix A as 


a = (n (m - Wd + i \...,n (m - Wd + i) ), 


with 


n {l) = (j> ({ 


kz, 


•d +1 


b l 


z d+ i£Z b i O ,ke{0,...,b l -1} 

Again analogously to |2j we obtain the following reduced fast CBC algorithm. 


Algorithm 2. a) Compute 0 (j?) for all r = 0,..., N — 1. 

b) Set r]i(k) = 1 + 710 ({^ 77 ^}) for k = 0,..., N - 1. 

c) Set zi = 1. Set d = 2 and recall that we have defined t = max {7 : Wj < m}. While d < min{s,f}, 

1. partition rj d _ 1 into b Wd vectors 77 ^,..., r 1 of length b m ~ Wd and let rj' = V^d-i + ■ ■ ■ + v!d-i * 
denote their sum, 

2. let T d (z) = , 

3. let z d = argmimTrfO), 

4 . let r] d (k) = %_! (l + 7 d 0 ({^r^})) for k = 0,... ,N — 1, 

5. increase d by 1. 


If s > t, then set zt = Zt +1 = ... = z s = 1. Then we have 


N -1 


R s n, 7 (b wl Zl ,b Ws z s ) = Ss(k). 


k =0 


Using pjj II2J HBJ CE] we find that Algorithm [2] has a construction cost of 

( min{s,t} 

A log N + min{s,t}lV + N ^ (m — Wd)b~ Wd 

d =1 


operations, in comparison to 0(sN log N ) operations for the standard CBC algorithm used for example 
in [ 20 ] . 


5 Conditions for strong polynomial tractability 

Let z = ( b Wl z \,..., b Ws z s ) € Z S N w be constructed with Algorithm|T]or[2]and consider the corresponding 
lattice rule. We are interested in conditions for tractability of the weighted star discrepancy of such 
lattice point sets. From m and we know 

uC[s] V V 77 


For now, let us assume that the 77 b w 'i ’s are summable, i.e. 
Sinescu in [0] and [20], we see that in this case ( 11 . 21 ) implies 


OO 

J2 7 jb Wj < 00 . Similar to Joe and 
3 = 1 


D* n „(z) < O 
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where the implied constant in the O-notation is independent of s. 
Recall that we have defined the information complexity as 

N*(e, s ) := min jiV € N : D* Nl (z) < e| . 


If 'ljb Wj < oo, it is easy to show that it is equivalent to consider the standard notions of tractability 
3 = 1 

with respect to N*(e, s ) or with respect to min jiV € N : R s Nry (z) < e|. 

Theorem 13.11 yields 


*Jr, 7 W < 4 n (ft + (l + 7,-ftv). 


3 =1 

We study the right-hand side of the latter inequality. 


4 n (%+(i+ 7i s„) < i p (ft+(i+ 26 —<»,.»•>) 7 . 2 (log 4 

3 =1 1=1 V V L 


+ 1 


< 


N 


i=i 


^ n (i + 7i (l + 4 (l + 26 min log iv)), 


where we have used 


j=i 


LfJ 


S »= E 


h^O 


h= 1 


/l 


IV 

T 


+ 2 < 4 log IV. 


n (ft- + (l + 26 min ^- m >) 7 j 4 log iv) (5.1) 


The second to last inequality is a well-known estimate for partial sums of the harmonic series. 
Now we have 


^ n (ft + (l + 26 min{u, j’ m} ) 7 ^) < 1 n (1 +7,' (1 + 4 (1 + 26^) log IV)) 
I=i I=i 

^na + ^iogiv). 

1=1 

oo 

Define < 7 ^ := 13 ]>4 7 jb Wj for d > 0. From [3J p. 222] or [61 Lemma 3] we know that 
j=d +1 

n (1 + 13 7j 6^ log IV) < (l + a^ 1 ) N^ 0+1 ^ d . 
i=i 


For 0 < <5 < 1 choose d large enough such that era < S +1 . Then 

S 

n (1 + 137 jb w > log AT) < c^ s N S , 

3 =1 

where c 7j( 5 is independent of s and N. Thus we have 

R s n ,j( z) < c^sN 6 - 1 . 

We obtain c-y^IV 5 ” 1 < e and thus R S N ^(z) < e if N > (c-y^e -1 )*^. Hence, if the jjb Wj ’s are summable 
we always achieve strong polynomial tractability. 
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Remark 2. Whether the conditions on the 7 j ’s and Wj’s can be mitigated while at least polynomial 
tractability still holds remains an unresolved problem. 
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